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Abstract 

We study induced nucleation by considering the accumulation rate of shrinking 
subcritical bubbles. We derive the probability for a collection of subcritical 
bubbles to form a critical bubble, and argue that this mechanism could well 
play a role in electroweak phase transitions if the Higgs is heavy. 
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Cosmological phase transitions have recently received much attention, in particular 
in the context of electroweak baryogenesis First order phase transition and bub- 
ble dynamics in the Standard Model have been studied in detail, and it has become 
increasingly clear ^ Q that for realistic Higgs masses, much heavier than 60 GeV, 
the electroweak phase transition appears to be only weakly first order. If the Higgs 
mass rriH ^100 GeV, both lattice studies and perturbative calculations run into tech- 
nical troubles. It is however conceivable that for such Higgs masses the electroweak 
transition is close to a second order one, and that it proceeds not by critical bubble 
formation but via thermal fluctuations. These would give rise to a considerable phase 
mixing at the critical temperature. 

In the context of the EW phase transition, subcritical bubbles were first discussed 
in [0. In it was shown that phase equilibrium can be reached provided the tran- 
sition is weak enough. In these papers, however, the disappearance of the subcritical 
bubbles was not accounted for. This happens in two ways: the subcritical bubbles, be- 
ing unstable configurations, tend to shrink; the bubbles are also a subject of constant 
thermal bombardment so that they may disappear simply because of thermal noise. 
The thermalization rate of small-amplitude configurations near the critical tempera- 
ture has been estimated in the EW theory in |^, where it was found that compared 
with typical first order transition times, thermalization is rather fast. Kinetics of sub- 
critical bubbles has been investigated by Gelmini and Gleiser [Q, who found, with 
a specific assumption about the form of the destruction rate due to thermal noise, 
that thermal noise becomes subdominant as the Higgs mass is increased. Their study 
deals with phase mixing above the critical temperature, and the possibility that pre- 
tansitional phenomena might play a role in the EW phase transition. Phase mixing 
has also recently been simulated numerically by Gleiser |^ in a 2+1 -dimensional 
model, with some interesting results. We shall return to the issue of bubble shrinking 
by termal noise at the end of the paper. 

In the case of a very weak first order phase transition subcritical bubbles may coa- 
lesce and form large regions of broken phase at and below the critical temperature, and 
effectively trigger the phase transition. This happens if the probability of subcritical 
fluctuation is large enough, in which case a collection of subcritical bubbles can form a 
region of the size of a critical bubble before they have time to shrink away. This might 
be called induced nucleation. Although it is clear that some dynamics is also involved, 
in the present paper we study the formation and growth of subcritical bubbles in a 
probabilistic approach. We derive an expression for the random growth of subcritical 
regions and compare the rate with the formation rate of critical bubbles. Interestingly 
enough, it turns out that for a large set of the parameter values induced nucleation 
is indeed possible, although for an extremely weak phase transition its completion is 
prevented because of the large size of the critical bubble. 



Conventionally two possibilities for critical bubble formation in d + 1 dimensions 
is considered. These correspond to the formation rates 
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Eq. (1) describes the situation where the critical bubble is formed by quantum tun- 
neling, and Eq. (2) by large over-the-barrier thermal fluctuations. Note the absence 
of h in Eq. (2). Sd[4>cr] is (i-dimensional euclidean action for the critical bubble con- 
figuration 0cr- Such critical bubbles are a result of a single, instantenous fluctuating 
event through or over the potential barrier. The timescale r^^ of such a critical bub- 
ble of radius Rc is determined by the thermal mass m{T) and wave number kc — ^/ Re- 



in a very weakly first order phase transition, however, the two phases are very 
much mixed already at the critical temperature, and thermal processes could create 
critical bubbles from small, unstable subcritical bubbles of broken phase. Instead of 
quantum tunneling or large fluctuations, the phase transition could be triggered by 
small semiclassical thermal fluctuations. This is called induced nucleation. 

Subcritical bubbles will of course disappear very rapidly, but if their production 
rate is large enough, the total volume in the broken phase may actually grow be- 
cause of the presence of fluctuating subcritical bubbles. Our aim is to determine the 
timescale in which a critical bubble is thermally built up from subcritical bubbles by 
calculating the probability for the appearance of a growing (spherical) subcritical bub- 
ble. Following we assume that thermal noise plays a subdominant part, and that 
the disappearence of the bubbles can be described in terms of shrinkage only. As we 
shall argue, this seems to be a reasonable approximation in our case. 

Let us assume that the subcritical bubble is a gaussian, spherically symmetric 
configuration with diameter /, i.e. 



where subcriticality implies that / ^ Rc- At finite temperature we need the d- 
dimensional action 



which essentially determines the formation rate of subcritical bubbles per unit volume: 



l/cu= [m{Tf + kl]-^'\ 
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Vv{l) = T[<Pi]/V c^T 



(5) 



In order to be spesific, we shall deal with a phenomenological potential for the 
order parameter (p suitable for a simple description of a first order phase transition. 



given by 

= ^m(r)V' - laT^' + (6) 

Most of the analysis takes place at the critical temperature (or slightly below), 
given by the condition 

miT^r = (7) 
At the critical temperature the non-zero minimum f (T^) of the potential reads 

v{T^) = IjT^ (8) 



and the correlation length is 



Ic = -JTfTT- (9) 



We shall let the dimension of the space-time to be a free parameter; the number of 
spatial dimensions is denoted by d. Note that this means that the parameters a and 
f3 are, in general, dimensional; only when d = 3 they are dimensionless. 

Given the phenomelogical potential Eq. (^, we can find out the action for the 
subcritical configuration Eq. (01) at Tc. The result is 
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The formation rate for subcritical bubbles is then given according to Eq. (^. 

In the cosmic soup there should exist all sizes of subcritical bubbles. We shall, 
however, use in our calculations a mean subcritical bubble I = Ic with an associated 
nuccleation rate Ty. With such an average subcritical bubble we are able to write down 
the probability for a process in which a large number of small subcritical bubbles are 
created adjacent to each other so that they form a critical bubble . 

Let us first find out how large a part of the space is occupied by the broken phase 
before the possible formation of critical bubbles, i.e. by subcritical bubbles in general. 
This problem has also been analyzed, in the case of critical bubbles, in []TU| and in [III 



however, without the shrinking effects (see also 0]). Let the radius of the subcritical 
bubble be Rsc and its life-time 

tsc = Rsc/V , (11) 

where v is the (average) speed of the shrinking bubble wall. Let us consider a volume 
V and, for the moment, discretize the time into infinitesimal intervalls of lenght St. 
We denote the volume occupied by the broken phase at the instance n6t (n G Z_|_) by 
Vb{n). Between the times {n — l)6t and n6t Nn = TvSt{V — Vb{n — 1)) new bubbles 
of broken phase have been fluctuated, each having the volume 7r'^/^-Rf^/r((i/2 + 1). 



Similarly, at {n — l)6t there were A^n-i bubbles but at t = n6t their volume has shrunk 
to Ti'^/^iRsc - v5tY/T{d/2 + 1). This procedure can be continued until all coexisting: 
bubbles have been counted. Correspondingly, in the broken phase N = TydiVbln — 1) 
bubbles of the symmetric phase are created, and so on. Note that at the critical 
temperature Ty = and V = V. Taking into account both effects, and taking the 
limit 6t 0, we obtain an integral equation for the Vj, to V ratio 

min{x, 1} minja::, p} 

^^s{x) = k f dy[l-s{x-y)][l-y]''-- f dy s{x - y)[l - y-]'' (12) 
V J P J P 

with the boundary condition s(0) = 0, where x is the scaled time x = tv/Rgc, p = 
{vR'J/{v'Rsc) and 

h — r 1"' — y-^sc) f^o\ 

^"^^r(rf/2 + i)t; ' " "^^r(rf/2 + iv • ^^^^ 

Here v' is the velocity of the unbroken phase bubble wall and R'^^ is the radius of 
the subcritical bubble of symmetric phase in broken phase. We have not taken into 
account the thermal noise effect described in It would modify the definitions of 
k and k' but here we assume that the resulting changes are small. From the integral 
equation Eq. (|12|) the asymptotic, equilibrium value of s can be easily solved: 

k 

lim s(x) = Sen = -, ; ;-• (14) 

^ ' ^ d+l + k + k' ^ ^ 

The equilibrium value at the critical temperature (where k = k') has the intuitively 
natural property that Seq \ when k oo (Fy oo). Equation Eq. (|T^) tells us 
how large a part of the spatial volume is occupied by the broken phase after the system 
has relaxed to equilibrium in an unorganized way so that no bubbles of critical size 
form. Note also that when the temperature decreases, k increases and k' decreases so 
that Seq increases, too. 

The issue at hand is then, when are the circumstances such that the subcritical 
bubbles can form a critical bubble? If this happens, the phase transition can be 
triggered by subcritical bubbles and completed by the expansion of the created critical 
bubble. Here we have to be slightly under the critical temperature because no bubble 
dynamics is really present exactly at the critical temperature. We assume that the 
critical bubble is built up layer by layer and consider spherical bubbles only, because 
they minimize the number of the subcritical bubbles needed for each layer and thus 
maximize the probability. Let us say that the radius of a bubble, still subcritical, 
which exists at the time t, is R. In order to that the bubble is really growing, in the 
time 2Rsc/v new subcritical bubbles have to fill a layer of volume Vr = 7i'^^'^[{R + 



2RscY ~ R'^]/T{d/2 + 1). Because the process is Poisson distributed, the probabihty 
that bubbles are nucleated is 

PN = ^A^e-^ (15) 

where 

A = TyVn^^ = 2k[{^ + 2Y - {^f]. (16) 

^ ^SC ^SC 

To fill the volume Vr at least Nr = VRT{d/2 + l)/{7T'^/^R'^J = {R/ Rsc + 2Y - {R/ R^cY 
bubbles are needed. The probabihty that the layer is filled is thus 

piR) = PN- (17) 

N>Nii 



A 

P{N, A) = — ^— / dtt^-^e-' (19) 



With a little algebraic acrobatics the formula Eq. ([TtD can be cast in the form 

p{R) = P{Nr, A) (18) 

where 

P(N W = 

r(iv)„ 

is a scaled incomplete F-function. 

We have derived the probability p{R) that a new layer of 'elementary' subcritical 
bubbles is formed fast enough so that the subcritical bubble size increases. We have 
neglected processes where large subcritical bubbles merge. This might be justifiable 
because a growing subcritical bubble should be a rare event so that the distance 
between any such configurations is large. In any case, such processes would only make 
the growth rate larger. 

When one uses a phenomenological potential of the type Eq. (^ there arises an 
extra complication because in realistic theories the order parameter (e.g. a Higgs field) 
often also has a phase. Then it is not enough that subcritical bubbles form next to each 
other but their phases should be correlated, too. Otherwise there should arise a domain 
wall between the two bubbles, and creating such a wall would require extra energy. 
In principle, one can easily account also for the phase by introducing a parameter 
fl G [0, 1] which determines the probability that the phase of the nucleated subcritical 
bubble is correlated with the phase of the pre-existing bubble. The probability for 
formation of N bubbles is then modified to read 

p'n = (20) 

and thus the probability of layer formation is 

p{R, n) = e-(i-^)"P(iVfi, nX). (21) 



The value of Q is, however, a dynamical question which we are unable to address 
here. When the new subcritical bubble overlaps with the pre-existing bubble, the 
latter may influence the phase of the newly forming bubble and may even force the 
phases to be correlated. In such a case Q = 1. Subcritical bubbles, which could be 
considered as wave packets of elementary quanta, might also have a velocity relative to 
each other. In that case energetics would no longer hinder the removal of the domain 
walls, resulting in fl = 1. In what follows, we shall write down our general expressions 
for an arbitrary fl, but when applying the results, we shall merely assume that Q = 1. 

We are now ready to write down the probability that a subcritical bubble is growing. 
It is simply the product over all layer probabilities, 

p{k) = l[p{Ri, n) = expE In P(i?,, n)] = expbX^, n, k)], (22) 

where i?j = {2i + l)i?sc and i runs from i = 1 to i = [1/2 + Rc/2Rsc] + 1, i.e. to the 
value where the bubble has reached the critical size. For large Rc/Rsc the function j 
can be approximated well by the integral 

1 

j ~ -JdRlnpiR^n) (23) 







Rs 



= ^ j du[lnP{N{u),2knN{u)) -2k{l-n)N{u)], (24) 



where N{u) = {u + 2)°^ — m*^. Now j can be computed numerically once the ratio 
Rc/Rsc, ^, d and k are given. 

To proceed further we have to calculate the average formation time r and the thick- 
ness (5 of a layer. Generally 6 < 2Rsc because the bubble is shrinking simultaneuously 
as the layer gets filled. N bubbles in the volume A^7r'^/^i?fp/r((i/2 + 1) are nucleated 
during the time t with the probability 

and at least N bubbles with probabihty 

Y: ^ (vkN±-) ^ e-^'^'t^ = PiN, vkN^) . (26) 

n>N ^sc'' ^sc 

Thus the average flUing time of a layer is 

t=oo 

{t)= j tdP{N,vkN^) = -^. (27) 



That is, the average time to form a new layer is Rsc/{vk). Because during the same 
time the bubble shrinks an amount Rsc/k, the bubble radius after z*^ layer can be 
solved recursively. The result is Ri = i{2 — l/k)Rsc- The critical bubble size has been 
reached when i = ic = {2 — l/k)^^Rc/Rsc and thus the formation time of critical 
bubble is 

. R^r Rr/V /„„s 

vk 2k-V ^ ' 

The average layer thickness is given by 

5 = Ri+i -Ri = {2- y)Rsc. (29) 

k 

To complete the calculation we have to write down the formula for the probability 
Pc that a subcritical bubble grows up to a critical one. Here we have to use the layer 
formation time r and thickness S as parameters, that is A = Fy (t) 7r'^/^[(i? + S)'^ — 
R'^]/r{d/2 + 1). The calculation results in the expression 

Pcik) = e=% (30) 

where 

' " l/kJo 

and Nc{u) = (n + 2 — l/k^ — u'^- Inspection shows that in the case f2 = 1 the value 



J- = o r"'^^ [lnP(iV,(M), nN,{u)) - (1 - n)N,{u)] (31) 



of jc, as a function of k, converges rapidly to the asymptotic function jc{Rc/Rsc) = 
jc{Rc/Rsc, ^ = I, k ^ oo). A plot of jc{Rc/Rsc) is presented in Figure 1. For Q = 1 a 
good fit for any A; > 1 is 

p / R \ 1-01+0. 008/(fc-0.2)i'20 

Jc(^,l,A:) = -10-°-^^-^°-^^A"^-^ (1^) . (32) 

For decreasing Q jc increases, however, rapidly even for as small ratios as Rc/Rsc = 10. 
In Figure 2 we demonstrate this behaviour by presenting function of for 

k ^ oo and Rc/Rsc = 10. 

We may now apply the rate Eq. (|34D to cosmology. Because Eq. (|34D states that 
Tyj oc e~^^^'^^'^'-^~^''\ we have to compare the combination pSd[(l)i] — jc = S^^^ with the 
Hubble rate. An induced critical bubble is obtained when 



M, 
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S^ff^\n[-^) . (33) 

The probability Pc{R) is, in general, very small. The number of subcritical bub- 
bles during the phase transition, i.e. the number of possible seeds, can, however, be 
very large. In this case the critical bubble formation rate is given by 



ry,/ =Pc{k) Fy. 



(34) 



where Ty can be computed by using Eq. (p!0|). To get some idea of the magnitudes 
of the parameters needed for induced nucleation, we compute the rate Eq. (|3^) in the 
thin wall approximation using the potential Eq. and assume that v is of the order 
of unity. 

First we need the size of the critical bubble. Assuming that there is only little 
supercooling, the bounce action can be written as @] 

^/^=)^^(X^-150, (35) 

where 

A ~ 1 - 0.0442—- = 1 - 6 . (36) 

Small supercooling, or 1 — A ^ 1 thus implies that a <^ 500 A^. Solving for A yields 
the transition temperature Tf. It then follows that at the transition temperature 

m\Tf) = '^-^{l-6)Tf . (37) 



Expanding the potential at the broken minimum (j) = v{T) = aT{l + y 1 — 8A/9)/2A 
we find 

- 6 = V{v, Tf) = -m\Tfy - -Xv' = -0.00218 Q + 0{6') . (38) 



The height of the barrier is situated at 0max — v/2 with V{(f)ma-Ki 



max 



a^T^^/(144A^). As T^. ^ Tf we may conclude that the thin wall approximation is 
valid if — e/\4iax = 0.314a^/^/A^/^ ^ 1, or a <^ lOA^, which is in accordance with our 
assumption of small supercooling. 

To get the size of the critical bubble we still need the surface tension. One easily 
finds 

Thus we finally obtain 

Rr 2a _ A3/2 



26.9^ . (40) 



Here Rgc — '^/f^iTf). Since a <^ lOA^ we see that Rc/Rsc S> 1 as it should. 

Combining Eq. ([10|), Eq. (^2]), Eq. ([33| ) and Eq. (|40|) at the critical temperature 
we obtain a condition for inducing a critical bubble: 

/ ,3s 0.506 

2.06-^ + 9.10 — <150. (41) 



A3/2 



a 



This implies a > 4.0 x 10^^ A^, that is, perhaps surprisingly, induced nucleation does 
not work for a too weak transition. The result can be interpreted so that in the range 
1 < A^/a < 250 where the phase transition is very weak, it is triggered by subcritical 
bubble formation immediately below the critical temperature. If the phase transition 
is weaker, the critical bubble radius becomes too large to be induced by subcritical 
bubble growth. On the other hand, whenever \^ /a < 1 the phase transition is not weak 
any longer: the radius of a critical bubble decreases and bubble dynamics becomes 
inportant so that our calculation loses its validity. The correlation problem connected 
to the parameter Q also becomes more important when the "seed bubble" is of same 
size as the newly produced ones. 

To appreciate the weakness of the transition required of induced nucleation, we 
may consider the 2- loop result for the electroweak effective potential calculated in Q]. 
Their result for = 87 GeV is well fitted by the potential of the type Eq. with 
a = 0.048 and A = 0.061. Thus in this case a ~ 210A^ is well outside of the induced 
nucleation range. As for larger Higgs masses a is likely to be smaller, it is conceivable 
that induced nucleation could play a role in EW phase transition. We should however 
emphasize that the result Eq. (^iD , and the subsequent estimates, should be considered 
as the most optimistic cases. Dynamics that correlates phases, and other issues related 
to the velocity of the subcritical bubbles or the surface tension between the bubbles, 
are likely to make induced nucleation more difficult, not facilitate it. 

Regarding thermal noise, we have compared the thermalization rate of small am- 
plitude configurations in the EW theory 0, given by F ~ 10~^T, with the formation 
time of an induced critical bubble, Eq. (^). We find that in the region where the thin 
wall approximation is valid, thermalization does not occur. This seems to indicate 
that thermal noise does not play an important role in induced nucleation. 

It would be of interest to compare the present approach with the results of Gleiser 

although it is not obvious to us how to do it. The problem appears to be a technical 
one only, though. Gleiser also employed a potential of the type Eq. (j^) (in d = 2) but 
simulated the heat bath by white noise which introduced an additional parameter to 
the problem. It is not clear how the parameters in his simulation are related to the 
parameters (including T^) here. 
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Figure captions 



Figure 1. The probability function j{r,Q = 1) defined in Eq. (|23|) in the asymptotic 
hmit A; — > oo as a function of the ratio of critical and subcritical bubble radii r = 

Figure 2. The function j{r, Q) in the asymptotic limit — > oo as a function of Q for 
the fixed ratio r = Rc/Rsc = 10. 
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